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1.  Introduction  and  Sunnary 

In  dealing  with  the  problem  of  the  reeponte  of  panels  under  blast 
loads,  such  as  from  the  firing  of  guns  or  rockets.  It  has  been  observe!  that 
the  severity  of  the  damage  depends  upon  the  speed  of  the  aircraft.  From  cer- 
tain elementary  considerations  of  the  nature  of  the  blast  load  (to  be  reported 
elsewhere) , the  reasons  behind  the  ptienomenon  can  be  qualitatively  explained. 

It  is  also  clear,  however,  that  the  surrounding  fluid  might  have  some  Influence 
on  the  response  through  its  reaction  on  the  transient  motion  of  the  panels  in 
question.  To  determine  the  approximate  nature  and  the  order  of  magnitude  of 
such  aerodynamic  effects,  an  idealised  two-dimensional  configuration  has  been 
studied  and  the  results  are  presented  in  this  report. 

The  simplified  configuration  treated  consists  of  the  following:  the 
skin  structure  is  represented  by  a flat  sheet  of  infinite  length,  over  which 
a uniform  stream  flows  at  a velocity,  \f  . The  motion  for  each  panel  is  con- 
sidered separately,  so  that  at  one  time  there  is  only  a "bump",  in  transient 
motion,  on  an  otherwise  flat  sheet.  The  aerodynamic  reaction  induced  by  the 
motion  of  the  bump  is  derived  on  the  basis  of  small  disturbances  and  incom- 
pressible potential  flow.  If  there  are  several  panels  excited  into  motion 
simultaneously,  the  resultant  aerodynamic  reaction  on  each  panel  can  be  con- 
structed by  means  of  superposition. 
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In  particular , a unlforn  panel  with  simply- supported  ends  la 
traatad  In  detail  in  terse  of  normal  coordinate*.  The  Important  dimension- 
t leee  paraaetere  entering  the  problem  are  the  ueual  reduced  frequency  and 
b&ii  ratio  (panel  to  air) , ae  in  all  oscillating  aeroelaetlc  phenomena.  The 
modifications  on  the  vibration  characteriatlce  are  discussed,  The  natural 
frequency  L each  eodc  tends  to  he  reduced,  ae  the  ar parent  case  and  the 
quasi- steady  forces  both  contribute  to  this  reduction.  The  reduction  in  the 
lower  natural  frequencies,  compared  with  the  frequencies  in  vacuo , is  of  the 
order  of  a few  percent  for  a typical  duralumin  panel  on  an  aircraft  soring 
at  an  average  epeed,  but  may  become  large  for  more  flexible  panels  and  for 
panels  on  submarines  under  water.  The  higher  natural  frequencies  are  practi- 
cally unaffected  by  the  surrounding  fluid. 

In  indication  is  further  made  for  the  condition  for  dynamic  inst*» 
blllty  based  upon  a representation  of  the  two  lowest  modes.  It  is  found 
that  for  typical  Aircraft  panels,  instability  is  essentially  caused  by  tho 
quasi-steady  forces.  The  aerodynamic  damping  makes  wary  little  contribur- 
tlon  under  ordinary  circumstances. 

It  Is  hoped  that,  despite  the  orer- simplification  of  the  two- 
dimensional  treatment,  application  to  actual  cases  may  be  made  for  first 
order  estimations.  This  may  be  done,  for  example,  by  using  the  concept  of 
the  "representative  section"  well-known  in  the  field  of  aeroelastlciiy. 
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2.  Method  of  Solution 
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Consider  the  complex  potential 

w(e).^e<ec,E 


(1) 


where  z * x + iy,  and  h may  he  complex. 

It  can  easily  he  shown  that,  within  the  scope  of  linearized  theory  for  small 
disturbances , such  a function  represents  the  disturbance  potential  due  to  a 
sinusoidal  boundary  of  period  , the  boundary  being  taken  to  be  close  to 

the  x-axls.  Regarding  Eq.  (1)  as  an  elementary  solution,  the  velocity  poten- 
tial for  any  arbitrary  boundary  shape  can  be  constructed  immediately  by  super- 
position, and  is  seen  to  be 


fo(ot)6  doc 


(2) 


or 


j'*4(ec)<Ll*'d0C) 


(3) 


Differentiating,  we  obtain 


<f) y ~^^{f  ^ e 1<>C  *dixr  J 


On  the  boundary,  therefore,  the  derivative  is 


$J  [*(-<*)  h (osje1**  d oc} 

J ti/  *0  L * 
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/: 


i**  T,.  .w-i** 


where  Jt(oC)  denote*  the  conjugate  of  hfl£) 


If  the  derivative  ^ / 

V ty*  O 


ie  expanded  hjr  the  Yourier  integral, 


we  nay  write 


K/U.t‘-  /"e  « “**  </< 


8f«)e  -,-Br-ftgJ 


-iw> 

® dot 


But,  with  the  relationship  that. 


6 I becomes 
M o 


y,.. x/ 


i<K  - -1**1 

+ d Coc ) e.  J clt 


Comparing  Bq.  (5h)  with  2q.  (ha),  we  may  conclude  that 

- oc h C oc)  * **  B(oc)  ’ (6; 

The  function,  B(w) , of  course  is  the  Inverse  transform  of  4<j/yo 

0| 

bw-17  / V«  dx  (?: 

-<0  • 

Hence,  with  the  given  normal  velocity,  • *h°  disturbance  potential 

can  he  evaluated. 
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3*  Application  to  the  Tranaient  Motion  of  a Bump 


In  studying  the  problem  of  the  transient  motion  of  a "bump" , the 
surfaces  adjacent  to  the  "bump*  are  assumed  to  be  flat  and  to  extend  to 
infinity.  The  air  stream,  of  Telocity,  1/  , is  taken  to  be  parallel  to  the 
flat  boundary.  The  x-axis  is  chosen  to  coincide  with  the  flat  boundary,  and 
the  "bump"  extends  0 rer  in  the  interval,  ^ 77  . Displacements  of  the 

■bump"  are  considered  small,  so  that  the  boundary  conditions  are  approximately 
satisfied  on  y 8 0. 

Thus,  if  y ■ f (x,t)  represents  the  Instantaneous  shape  of  the 
bmp,  the  boundary  condition  to  be  satisfied  by  the  disturbance  potential  is: 


Also,  in  order  for  the  disturbance  to  die  out  at  large  distances  from  the  panel. 


the  conditions 


< P and  (f)^  — 0 as  / x/  or  <j>  or  both 


■co  (9) 


is.. 


dx  'O 


must  be  satisfied.  Condition  (9)  1*  self-evident.  Condition  (10)  represents 
the  "closure"  requirement,  i.e.,  the  instantaneous  streamline  must  adhere  to 
the  flat  boundary  after  passing  the  bump. 

It  is  obvious  that,  if  the  instantaneous  shape  of  the  "bump"  is 


defined  by  the  equations 
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• f (*  **  ) 


Ol  X «Tf 


elsewhere 


The  shape  may  he  redefined  In  terms  of  a Fourier  series,  so  that 


(/  « »n(-f)  Sin  hk 


o s xSf 


elsewhere 


Here  It  Is  assumed  that  the  hump  starts  from  sero  at  x * 0 and  x ■ ir  . If  the 
"hump"  represents  a panel  of  homogeneous  material  and  of  constant  thickness, 
each  Fourier  component  corresponds  to  a normal  mode,  and  q^t)  corresponds  to 
the  transient  deflection  in  terms  of  the  nonnAl  mode  components.  By  introduc- 
ing Eqe.  (11)  and  (12)  into  Bq.  (8)  the  boundary  condition  becomes 

Vt.o’V  VnC**  Sin  °* 

= £ An(t)  cosCnw-^  , O 4.  x (13) 


H-p  ’ 0 


elsewhere 


And  consequently,  the  problem  is  reduced  to  finding  a velocity  potential 


satisfying  each  term  of  Bq.  (13),  namely, 


^/(^-O  C^n(f ) c°s(r7X-SnV) 


04  X 


elsewhere 


57 r) 

iere  J 


U3*> 


together  with  conditions  (9)  and  (10).  The  complete  solution  is  then  ob- 
tained by  superposition. 


Proceeding  as  indicated  above,  for  Bq.  (13a)  we  obtain  by  using 
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a (5)  -/J  ^l± — 
Sl (5^* 5th  _ — 


Both  integrals  in  Sq.  (16)  show  singularities  within  the  range  of  integra- 
tion, hat  difficulty  will  only  he  encountered  at  the  origin,  of  * 0. 
Expanding  Xq,  (16)  into  partial  fractions  and  integrating,  ve  obtain 


tn/uno*  -^fcoS  r-  CM 


d±n  s in  /n[”c.os  n ^ I* 

Jfi  si£  <fa  £iin  j-cosnir  Ci((Tr-x)€)*Ci(*f  )} 
TT  n f-o  [ J 


Tor  even  values  of  n,  ve  have  (Ref.  2) 

fiZ  [-  C ("•*)()+  O JL/gz  / (20) 

Bat  for  odd  values  of  n,  there  is  a divergence  of  <pn  even  for  large  values 

of  j t(  , and  thus  condition  (9)  is  not  satisfied. 

The  source  of  this  difficulty  can  he  tr^tced  tc  the  fact  that  the 

closure  condition,  condition  (10),  is  not  satisfied  by  that  part  of  4qL%0 

® d 

specified  by  Xq.  (13a) , which  is  associated  with  sin  /n  . 
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Defining  the  function 

6 L,0~An  Sen  n*  Sen  /„  y o ± x * TT 

/ J 


elsewhere 


we  eee  that 


r"  • r 

J Aft  Sin  fn5<-n  nx  d*  (/-cos  n'Tr) 


\ o if  ri 


la  odd 


(21) 


Zq.  (21)  definee  a constant  displacement  of  the  stream  line  In  the  y-direction. 
As  shown  In  Appendix  A,  the  addition  of  another  "hump"  in  the  form  of  a step 
eliminates  this  constant  displacement,  and  the  closure  condition  is  satisfied* 
The  dlrergent  term  in  Zq.  (19)  then  disappears.  It  follows,  therefore,  that 


where 


d)  I a ~ An  cos  nlT  J/fr-x,  "J-X,  (xj  r>)l+ 

~TT 

^ n Si  n ^["cos  hir  (7T-X,n)-X2(x,n)+IJJ 

(19a) 


^3=  ^ for  even  t) 

- 0 for  odd  r> 


(22) 


Rewriting  I,  and  Ig  according  to  their  definitions,  Zqs.  (17).  we 


obtain  finally, 
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K / ‘ n * ^ <r-n x)  c"  nx]* 

*JLl[c,(nx)  cos  nx-S,  fnx)  5/*  »xj+ 


l + Cts  nJJT 

Tv  n 


A„s  ^J^rJ 


(23) 


where  Anc  and  are  abbreviations  for  cos  f and  sin  6n 
and 

C |(nx)'-  Ci  (nV- nx)  - Ci  (r?  x) 

S,  (r»x)  ~ St  (wTT-nx)  + Si.  ("*) 


re spectively , 
(24) 


Since,  by  definition, 

c<  (-})-- Ci( -3) 


Sif-p- 

the  condition  that 

(A»^u«o  0 ft9  x~~“’ 

is  satisfied  as  required. 
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i*.  The  Pressure  Distribution 


Baring  obtained  the  solution  for  f a in  Iq.  (23) , vs  may  now 


evaluate  the  pressure  distribution  along  the  "bung”.  The  pressure  is  given 


» - e u ~zk/ 

A ' T • ^ tO 


with  the  pressure  on  the  underside  (7  3 -0) , taken  to  be  the  same  as  that  at 


infinity. 


Differentiating  Eq.  (23).  we  obtain  for  the  n"1  component. 


'V 


— 


“jfC  x)J sin  r?x  + 

rS'(nx)4.Ct(n>)'lcoi  r,x)  + 


r 

^nxifc!Cnx)-S,(nx)]co5  n* - 

IT  ( L 2 

C»*;7S<n  * 


1+  Co s T>  IT  ^ _ _/ j) 

ZYT n 1,5  t /*/  /V-x/  j 


where 


' J C,  x ) - - k] 


(n ->  - ^ f ^ ^ J* 


After  simplification,  in  the  region  OiX-'TT  , i.e. , along  the  "bump* 
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-US.  Jc,  Cnx)  cos  nv  — S/ Cnx)  s/#7  n *]  — 

^ * 


^11.]"  c,  (nx)  Sin  vx  -t*  ^(nxj  c*s  ox  -+  (28) 


l-C0S0'TT/r»  f \ -7 

-XS )] 


Tor  points  beyond  the  "hump”,  a corresponding  expression  may  readily  he 


written,  which  should  he  used  when  the  mutual  influence  of  two  adjacent 


"humps1'  is  to  he  computed.  This,  however,  is  omitted  here. 


Since,  hy  definition,  the  following  symmetry  holds: 


C,  (nx)-*  - C , C O 'Tf-n  *) 


Si  (n*)«  S,(nlT~rix) 


it  may  then  he  verified  that: 


(1)  for  n * odd,  the  part  of  $njXJ  . Q in  Iq.  (28)  due  to  Anc  is 


symmetrical  with  respect  to  the  mid-point;  the  part  due  to  AQg  is 


anti- symmetrical; 


(2)  for  n • even,  the  part  due  to  AJ1C  is  anti- symmetrical  with  respect 


to  the  mid-point;  the  part  due  to  Afia  is  symmetrical. 


It  will  he  seen  that  the  last  term  in  Sq.  (28)  involves  singulari- 


ties at  the  leading  and  trailing  edges.  This  may  he  interpreted  as  being  due 


to  an  effective  movement  of  the  segment,  CSx<  'TT  , as  a rigid  body  from 


the  flat  boundary.  The  rest  of  the  terms  result  from  the  deformation  about 
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this  means  position. 

For  the  contribution  due  to  ^ , the  same  expression,  Iq,.  (23)  • 

may  he  used  except  that  the  amplitudes  A^  and  Aq>  must  he  replaced  by  their 

derivatives  A^  ■ -flt  (An0)  and  kn%  * ct_  (An#)  respectively.  Thus 

dt  dt 

<Pn±l  * — 222  fC,  (ox)  S in  nx  + 5,(r»t)  cos  njeJ  + 

J nir  L 


Aos  fc^nx)  CoS  nx  HX  + 

orr  *- 

<»> 


I -f-CQS  n TT  Q_  / x 
2. 


Ve  are  now  in  a position  to  evaluate  the  total  change  of  pressure 
due  to  the  motion  of  the  "bump " t 


i *./  -fi  <> » / 


with 


*mL  “4  f„,t  4 


given  by  Sqs.  (28)  and  (29). 


(30) 


5.  Generalised  Forces  on  a Uniform  3 imply- Supported  Panel  IXie  to  Aerodynamic 
Pressure 

For  a homogeneous  panel  of  uniform  thickness,  the  Fourier  components 
of  the  deflection  are  identical  with  the  normal  modes.  Tho  coefficients, 

(^(t) , in  2q.  (12),  therefore,  represent  the  deflection  in  terms  of  normal 
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coordinate*.  We  may  then  write  the  equations  of  motion  ae 

ir  “V  (3l) 

where 

Kj.  * generalised  mass  for  the  r**1  mode, 
qp  • deflection  of  the  r**1  mode, 
uj  s r*11  natural  frequency, 

Qp  - generalised  force  exciting  motions  in  the  r^  mode. 

It  Is  well-known  that  Qp  may  he  derived  from  the  principle  of  virtual  work. 
When  a pressure  distribution  acts  on  the  panel,  the  virtual  work,  &W,  given 
by  a virtual  displacement,  6«^  , is 

&W  x f - A p die  S ^ • (32) 

Impressing  y in  terms  of  normal  coordinates,  we  find  that 

8’S,  S"Si"  "* 

Sir*  r?x 

0 n*i 


Substituting  into  Xq.  (32) , we  obtain 

sn 


/n 

-Ap  sin  r*  < ix 

«!*•  0 


(33) 
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The  ten  My  la  glean  try 

Mr*  J*- j S loarxdx 


where 


yM  a naae  per  unit  apan 
N s total  naea  of  the  panel 


With  Iq.  (30)  for  A p,  It  follova  that 


• 


<? f f>  „/  Sin  r*dx 
ft*/  4 rnix/«j*o 

+ ^I  f A / Sin  rjcix. 
u*i  Jo  nitJ^%0 


In  ualng  Iq.  (35)  we  mat  take 


^h5  "*  r»  > A 


lnc*  Vn<\n 


aa  defined  h j Iqa.  (13)  and  (23).  Tor  abbreviation,  Iq.  (35)  nay  now  be 
rewritten  In  the  form 


-L 


- S.  -IT  Ji-WL  <J  „ 5 IniT.  + ? T *"<•-  t ' 

ft*»  or*  v o »»  rr  v 


and  o(  , ^ are  identified  aa  followat 

"in  rjf[c,c».)ces  nx-S,(ox)sm  nxf  t+SZLIlJLXn^^dx 

«[*>*>  S*"  ',^S'fn<)COS 

» - ££  r.w^C|0>*)ces  nx-S,  Cr*x.)Stn  nxjd*  J 


AERO-ELASTIC  AND  STRUCTURES  RESEARCH 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
Department  of  Aeronautical  Engineering 


PAGE  16 


CONTRACT  NO.  15  ori-07833 
HR  064-259 


Fhyelcally , 0 repreeente  the  apparent  maaa  effect  of  the  node  on  the 
r**1  node,  ^ n,r  repreeente  the  correeponding  damping,  and  y repreaentt  the 
quae i- eta ady  aerodynamic  force. 

Ve  nay  now  Introduce  the  functlone, 

Js  (n,m)*  T^S,  (n>)  c°s  excise  l 

[ (38) 

TC(  (r),m)-£'rcl(rtx)si'>  n>x  d* 

Jc  f C,(t»x)c*s  rnxdx  J 


which  are  eraluated  In  Appendix  B. 

Baring  defined  the  functlone  of  Jq..  (38) , ve  may  write  Iq.  (36)  ae 


r~  n)  +£  + 


11^2/1.,  „ ..  tx 
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(n>r- n )-  * r-rr*)  + 

|-TS|(n,r-n)4-£ 

r«(i-i  A dx7 

n J V*  TT-X  1 J 

0 

tfn.r  = " ^rf 7 Jc/  r'  ">  (». r+  n>  ~ 

jT5j(n,r-n)+iJ^(t.lr*n)J 


Making  use  of  formulae  (B.l)  to  (B.12),  we  may  write,  for  n^r 

uCn>r’  £fr‘l~?r^t  [ l+  c‘5  <■•■<■  n>ir]  p-‘  (rtr)-C;  (nir)*  in^-J+ 


- — (-  urCM  r'T-)[^  '*?-C'(rirj]J 


where  ^ * Euler's  constant  * 1.781072. 


Similarly,  for  n * r,  o(n  becomes 

jpf'-n  Si "*-•)-  n-(^n 

Jj  (<il) 

In  a like  manner,  ^n,r  W b®  written,  for 

* » Si  3t 

(L  (42)  ' 
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and  it  ia  eeen  that,  for  n • r,  §nn*0. 

finally,  the  expression  for  Vn<r  may  he  written,  for  nfr 

y„,r*  - 3r(-  &p)(>  ♦ n)ir^[Ci  M)-Ci 

and,  for  n * r,  £n>n  becomes 


Va  now  have  all  the  coefficienta, e(  , $ , ^ 'a,  for  evaluating  the  generalized 
force,  Qj.,  aa  defined  hy  Eq.  (36).  It  should  he  noted  that: 

(1)  ®<n>r  v*o  only  when  n and  r are  both  even  or  both  odd,  i.e. , the  odd 

modes  do  not  contribute  to  the  apparent  rnaao  of  the  even  modea  and 
vice  versa; 

(2)  $ri,r  ^0  only  when  n and  r are  not  both  even  or  both  add,  i.e.  the 

odd  modea  contribute  to  the  damping  of  only  the  even  modee  and  the 
even  modea  contribute  only  to  the  odd  modes; 

(3)  ^>n  o c 0 i.e. , there  ia  no  damping  of  the  mode  due  to  its  own 

motion; 

(h)  Yr  n s ¥h,r  *•••  <iaa*i-»teady  forces  are  reciprocal; 

(5)  |(r  0 only  vhen  n and  r are  both  ever  or  odd,  i.e.  the  odd  modes 

do  not  contribute  through  their  quasi-steady  forces  to  the  motion  of 
the  even  modes  and  vice  versa; 

(6)  The  mutual  effects  of  neighboring  modes  are  in  general  larger  than 
those  of  widely  separated  mode?,  owing  to  the  factor,  r*-n*  , 
appearing  in  the  denominators  of  each  coefficient; 
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(7)  Tov  a glean  r,  the  apparent-mass  and  quasi-steady  force*  due  to 

j Qfj  fj 

other  modes  diminishes  in  the  manner— — for  large  n » r,  and  as 
0 r 

— - — for  large  r » n. 


6.  Iffects  on  the  Ylbratlon  Character! e tics  Doe  to  the  Aerodynamic  Pressure 

Combining  Eqs.  (31)  and  (36),  we  now  have  a system  of  an  Infinite 
number  of  linear  differential  equations: 


Lr~^r*»r)  S»r  ~ 


00 


0 

(WO 


where 


& n r = Xronecker ' s symbol 


~ 1 if  n * r 


= 0 if  n + r 

, . M 

K = ^2  = mass  ratio 

A new  set  of  natural  frequencies  may  now  be  found.  Assuming  that 
in  Hq.  (WO , 


- (cu't 

pm<t*e 

where  Qq  and  o'  are  the  magnitude  and  frequency  of  the  motion  in  the  n 
noroal  coordinate,  it  follows  that 


th 


4 
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r- 1,2 


The  new  natural  frequencies,  , therefore,  satisfy  the  Infinite  determinant 

/Afirf^O  '»»""«  (46) 

where  \ 

A _ * - /l'*(X  /„r—*Cir  )~‘  ?»r  tar”  Yhr  ) 


. y (“6«) 

A = modified  reduced  natural  frequency 
= TJ 

A.  = original  reduced  natural  frequency  ' 

U 

In  evaluating  the  infinite  determinant  it  is  necessary  to  consider  a finite 
number  of  modes.  It  is  obvious  that  since  the  aerodynamic  forces  are  smoothly 
distributed,  their  effect  is  likely  to  be  confined  only  to  the  first  few 
modes.  Tor  example,  let  us  examine  only  the  three  lowest  modes. 


-n\,3-X,3 


-iJl'Q. 


-si  YK-KjJ-x-a’-y.*  -iJi't 


Jl'  aC  ~ r 

31  9 31 
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It  la  easily  seen  that  the  result  will  be  an  equation  Involving  real  numbers 
only,  since  the  real  and  imaginary  terms  always  occur  alternately  in  the 
determinant.  For^j  * 0 a cubic  in  JX  is  obtained,  yielding  in  general  the 
three  modified  natural  frequencies.  The  solution  is  a function  of  two 
parameters  — the  mass  ratio,  K » and  the  reduced  frequency,  JX,. 

From  the  behavior  of  the  coupling  factors  between  the  n™1  and  the 
r modes,  it  can  be  shown  that  the  contribution  due  to  the  very  high  modes 

Q n 

is  at  most  of  the  order,  — — . Thus  for  the  lowei^-node  natural  frequencies,  the 

very  high  modes  probably  may  be  dropped  without  serious  error.  $uch  of  course 

has  always  been  the  intuitive  approach  for  systems  involving  a large  number 

of  degrees  of  freedom.  Assuming  the  validity  of  this  hypotheeie,  we  can  make 

the  further  deduction  that  the  very  high  natural  frequencies  may  also  be 

regarded  as  independent  of  the  lower  modes.  The  argument  is  as  follows. 

Vhen  the  r**1  mode  equation  is  written  down  for  r »1,  it  is  seen  that  the 

very  low  modes  have  contributions  of  the  order,  ^ • Since  the  very  low  modes 

n 

have  been  solved  from  a sub-determinant  involving  an  adequate  number  of 
modes,  their  appearance  in  the  equations  for  the  very  high  modes  may  be 
treated  as  being  an  external  forcing  function,  and  should  not  influence  the 
natural  frequencies.  Consequently  only  an  adequate  sub- determinant  involving 
modes  of  the  same  order  as  the  rtb  mode  is  needed  to  evaluate  the  rth  natural 
frequency.  Besides  “V  , only  ^ s exist  for  large  r,  and  are  of  the 

op4‘r-  i^l  • 

In  fact,  the  diagonal  terms,  are  the  dominant  terms  under 

ordinary  circumstances,  when  the  aerodynamic  reactions  are  relatively  small. 
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Without  the  aerodynamic  reactions,  all  terms  except  those  on  the  diagonal 
vanish.  Nov  for  large  values  of  r, 


A,r*=-HJll2-1.HJlr-Yrr 

Since  Jl*»’*of  the  order  of  r*  » «nd  Yrr  is  of  the  order  of  r,  while  $nr  is 
at  most  of  the  order  unity,  the  constant  parts  in  the  diagonal  terms  become 
very  large  in  comparison  with  those  in  terms  off  the  diagonal.  It  is  apparent 
that  if  SI*  ie  evaluated  from  the  diagonal  terms  only,  an  adequate  approxima- 
tion for  the  higher  modes  may  be  obtained. 

Th*  correction  to  this  approximation,  if  desired,  can  then  be 

* 

worked  out  by  a procedure  developed  by  Lord  Rayleigh:  Assuming  that  the 

system  is  excited  in  the  modified  r**1  mode  (which  is  almost  the  original  rth 
mode,  since  the  coupling  terms  are  assumed  to  be  small),  and  neglecting 
second  order  terms,  we  have: 
from  the  equation  for  the  n^  mode. 


^On  + Arn  ““  0 


from  the  equation  for  the  r^11  mode, 


T Anr  <fn  -f  Atr  fr  *0  W 

nxi  ® 

/ 

where  indicates  that  n = r is  excluded.  It  follows  inmediately , then, 

Arr  ~ £ Am  (lt, 

A 

Hnh 


* Ref.  3,  pp.  113  - 115.  also  pp.  136  - 137 
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The  right-hand  aide  of  Xq..  (^9)  is  already  a eecond-order  quantity.  Tor  a 
first  approximation , therefore, 

Arr*-(*-<{rr)-ri  +K  0 

whence  we  eee  that 

ji,2r  (50) 

r'  K“Xrr  r!  H r ‘ 

The  value  of  JV^  can  then  replace  the  unknown  fi*  in  the  right-hand  aide  of 
Xq.  (**9)  • and  a second  approximation  la  obtained. 

f ' (iC«„ -ui'..  g.  -r„, 

i IR,r  'i  (M. )(*.«„ 

(50a) 

Thus  the  correction  terns  are  mainly  those  due  to  the  neighboring  modes  n^s-r 
Tor  r7?l,  the  important  result  here  la  that  the  modified  natural 
frequency  is  given  by 

u)'r 


u)r  ‘ 2K-nr  ) 


(50b) 


81nce  y is  always  positive,  the  natural  frequency  is  reduced.  The  percen- 
tage reduction  decreases  as  r“^,  and  may  be  neglected  for  the  high  modes. 

The  reduction  Is  also  inversely  proportional  to  the  mass  ratio  of  the  panel  to 
air.  As  altitude  increases,  K increases,  and  the  reduction  in  the  natural 
frequency  becomes  smaller  and  smallsr. 

It  may  be  noted  that  Eq.  (^9)  and  the  succeeding  approximation 
procedure  say  also  be  used  for  the  lower  modes,  provided  that  the  diagonal 
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terms  remain  dominant.  Such  would  be  tho  case  when  ^ and/or  JV*,  is  large 
In  comparison  with  the  , l^nr  and  ^nr  coefficients. 


I 


7.  Further  Investigations  and  Some  numerical  Results  Based  on  the  First  Two 
Modes 

It  Is  desirable  to  be  able  to  Indicate  the  order  of  magnitude  of  the 
aerodynamic  effect  on  the  vibration  characteristics  of  the  lover  modes, 
following  the  general  remarks  on  the  high  modes  in  the  previous  section.  For 
this  purpose  let  us  consider  only  the  first  two  modes.  The  coefficients  of,  , 
etc.  are  found  to  be 

> <2  '-Kj,,*6  » 

« O , -/.»«,  $2)-/38£  , ~o  > 

Thus  the  frequency  equation  becomes 

£.JI'  ( -BA3][-  -A'  (k  + .22Z  )+kJl\ -S.6&] 

+ JL'Z(-l.<jSJl)(l.58£)  -0  (51) 

Noting  that  for  simply- supported  ends, -f-il ( , we  obtain  after  algebraic 
manipulation. 
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-fl'4 (Kf .£4i) (K  + .222 )-  f 17 H + 4/6  ) K Jl%  8.»f  « * 0.7*  Jj 

+ l(,(kSl*-.$S*)(KJL*-2+$)”0  (51a) 

Heglecting  aerodynamic  forces  correspond*  to  the  case  of  )(•»» , for  which  the 
two  roots  of  JT*  reduce  to  JL*  and  16J1*  , as  required.  Tor  large  values  of 
K , the  roots  can  be  expressed  In  descending  powers  of  K , in  the  manner 


Wi'j  and  being  the  modified  first  and  second  natural  frequencies  resulting 
from  the  presence  of  the  surrounding  fluid.  It  should  be  pointed  out  that, 
even  for  the  first  two  modes,  the  natural  frequencies  still  tend  to  he  reduced. 
Tor  the  higher  modes  the  same  conclusion  was  reached  above. 

When  JT*  is  large,  2q.  (51a)  may  be  simplified  by  dropping  terms  of 
the  order  of  • ^as,  for  fi*  ,»/ 

(Ji-)* (K+jt**XK *■***)-($■)  *(*“  ♦'**)♦»* *** o <»> 

/ i 

. i 

Tor  a first  approximation,  the  frequency  ratio,  , is,  then,  a function 
only  of  the  mass  ratio,  V . The  physical  interpretation  of  this  is  that  the 
modification  is  now  mainly  due  to  the  apparent  mass  of  the  fluid  associated 
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with  the  panel  vibration.  In  contrast , the  previous  correction  for  the  case, 
K >>  1,  bajt  be  regarded  to  originate  largely  from  the  quasi-steady  aerodyna- 
mic forces. 

Tor  numerical  examples,  let  us  consider  some  typical  values. 

(1)  Duralumin  panel,  12*  between  supports,  .064"  thick,  sea  level 
atmosphere , U a 50O  ft/ sec. 

jt.  r*  .026  slugs/ft 
G),  sm  870  rad/sec 
X e~  10.9 

Jl 

500  * 


From  Sqs.  (52),  we  obtain 


' vd 


(£> 

u), 


0.90 


( 0-97 

to* 


l.e. , 


^ — JS-  , 

(2)  8tesl  panel,  2^"  between  supports,  thick,  in  water,  U - 30  ft/sec 


•i*  — 39 

MJ, 


32  elugs/ft 
V*i 2,600  rad/sec 

X .08 


a **  2fe00S2  e 180 

* 3o 
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Xq.  (53)  "becomes 

A ot  f ■£•')* *.+*  ( l.ee  *« 

from  which  it  follows  that 

OIL  *»+  So  , iff  ***  .S’d 

to,  vPz 

A remark  may  he  inserted  here  regarding  the  applicability  of 
formulae  derived  in  this  report  to  panels  in  water  and  moving  at  relatively 
alow  epaede.  The  theory,  of  course,  is  still  valid,  provided  only  that  the 
deflections  he  sufficiently  small  for  the  linearization  of  the  boundary 
conditions.  The  typical  values  in  the  above  example,  on  the  other  hand, 
seem  to  indicate  that  \(  will  he  small  and  large.  If  K is  small,  the 
natural  modes  in  vacuo  are  poor  approximations  to  the  true  mode  shapes. 

Strong  coupling  effects,  therefore  must  he  expected.  Probably  more  than 
two  modes  are  necessary  even  for  a rough  estimate  of  the  modified  first 
natural  frequency. 

Some  discussions  on  the  possibility  of  dynamic  instability 

of  the  panel  may  also  be  made  on  the  basis  of  the  two-mode  representation. 

Let  us  consider  again  Xq.  (51a).  It  is  apparent  that  a divergent  motion 

£ 

sets  in,  if  the  value  of  Jl/  satisfying  Zq.  (51a)  becomes  negative  or  complex. 
A locus  in  the  K-  ilf"  , plane  can  be  drawn  to  separate  the  stable  and  unstable 
regions  in  terms  of  these  two  parameters.  Regarding  Xq.  (51a)  as  of  the 
typical  ferm 

A tC  = O 
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vr**  #**  t 


r?-5 


Ve  Me  that  the  condition#  for  positive,  real  roots  of  x are 

(1)  a,  r,  c >o 

(2)  B2  ^ 4AC 

8ince  K>9,  the  coefficient  of  J1  is  alvays  positive.  Hence  the  conditions 
for  stability  reduce  to 


(5*0 


(55) 


(a)  RJX*  2 8IIK--7* 

1 >?  K + +U 

£ 0.4?7*0(^-)  f.rK>?l 

(b)  ( K-a*  ~ 3 55  )(*-#.*- 2. 4*)  2-0 
i.e.  ,355*  £ K-fl?  or  KIL*  * £ +2 

<c)  fO?****.!*)*-^-8  % 

4*(k*.a4C)0<-K222  )(KA%2>3)(k-A*-  35-5) 

(5«) 

Instead  of  tracing  the  complete  boundary  defined  by  El.  (56) , a simplified 
criterion  nay  be  derived  for  the  cae*  El-  (58)  is  then  approximated 

TO  .*  77+  'Iff 

Within  el  Ids- ml.  Accuracy,  the  result  Is  (5^*) 


(57) 


Combining  the  conditions  (54),  (55)  and  (57) : we  conclude  that  for  U »l  . 
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stability  prevails  if 

XJl*  >2.43  (58) 

Physically,  this  means  that  the  quasi-steady  forces  tend  to  reduce  the 
effective  spring  constant  in  each  mode.  The  strongest  influence  is  on  the 
first  mode.  When  the  effective  spring  constant  of  the  first  mode  is  reduced 
to  a negative  value,  divergent  motion  will  naturally  result.  The  same 
influence  on  the  higher  modes  is  much  smaller.  ForK>:>J  * furthermore,  the 
damping  and  apparent  mass  terms  are  of  no  consequence  in  the  determination 
of  the  stability  range. 

For  a numerical  example,  ?„et  us  take  a duralumin  panel  of  .064-H 
thick  and  46n  between  supports,  at  seat-level  atmosphere.  One  finds 

/k.  ,026  slugs/ft 

U)(  55  rad/sec 

K 2.7 

2 

The  critical  value  of  JT,  is,  by  lq.  (58) 


r * 

Cl*.  2.7 


Hence,  the  "criticeil  speed11  is 

(-ft- 1 'c  r. 


%C. 


It  must,  however,  be  stressed  that  the  results  in  this  section  can 
be  valid  only  when  a two-mode  representation  is  adequate  to  describe  panel 
vibration.  The  discussion  on  stability  is  further  restricted  by  the  assump- 
tion of  small  disturbances.  Whether  the  stability  involving  finite  deflec- 
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cions  haw  roughly  the  seme  criterion  or  not,  ie  a question  beyond  the  scope 
of  linearised  theory.  It  may  further  be  noted  that  dynamic  instability 
described  abore  is  different  from  the  conventional  "flutter”  phenomenon.  In 
■flutter",  ws  may  recall  that  the  resultant  damping  of  the  system  vanishes 
at  the  critical  speed,  while  here  the  resultant  spring  constant  is  the 
criterion,  at  least  for  K » I (relatively  heavy  panel). 


8.  Summary  of  Results. 

(1)  Based  upon  the  linearized  theory  for  small  disturbances,  the 
velocity  potential  of  arbitrary  transient  motions  of  a two-dimensional  bump, 
on  an  otherwise  flat  surface,  in  a uniform  incompressible  stream  has  been 
determined  in  terms  of  the  Si  - and  Cl  - functions. 

(2)  Application  of  the  theory  is  made  tc  the  problem  of  a simply- 

supported  uniform  panel.  The  effect  on  the  vibration  characteristics  depends 
on  the  usual  parameters,  (the  reduced  frequency  of  the  first 

natural  mode)  and  (the  mass  ratio  of  panel  to  air).  Coefficients 

representing  the  apparent  mass,  damping  and  quasi-steady  forces  are  obtained, 
in  a form  for  use  with  normal  coordinates,  which  are  strictly  true  for  a 
simply- supported  panel  in  vacuo. 

(3)  Tot  the  higher  modes,  only  the  quasi-steady  aerodynamic  forces 
have  some  contributions.  The  result  is  a lower  natural  frequency.  The 
percentage  reduction,  however,  decreases  very  rapidly  for  higher  and  higher 
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modes.  It  Is,  furthermore,  inversely  proportional  to  the  mass  ratio,  X , 
and  the  square  of  the  reduced  frequency , -0-f. 

(4)  An  iterative  procedure,  following  Rayleigh,  is  Indicated  for  the 

evaluation  of  modified  natural  frequencies  whoa  X and/or  . Panel* 

on  aircraft  usually  belong  to  the  category  of  large  K . while  panels  on 

r\2 

submarines  in  water  may  correspond  to  the  latter  case  of  1LI  »|  . 

✓ 

(5)  Panel  behavior  based  on  a representation  of  the  first  two  modes  is 
studied  in  mors  detail.  Approximate  formulas  are  given  for  cases  when  K»l 
and  when  >> ! . It  is  again  found  that  the  natural  frequencies  tend  to 
diminish  in  the  presence  of  the  surrounding  fluid. 

(6)  A criterion  for  dynamic  stability  is  also  derived  based  on  the 
two-mode  representation,  namely,  for  stability 

Ufl*  >2.43 

To  illustrate,  a duralumin  panel  of  4*  span  and  .064*  thick  is  estimated  to 
have  a critical  speed  of  about  230  ft/sec  at  ses^- level. 
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APPENDIX  A 


Modification  on  for  Odd  ft  , to  Satisfy  Conditions  at  Infinity 


Let  us  consider  the  hump 


0^/  *AnSirt”x  sin  for  o<x<'7Tl 


(Al) 


- q elsewhere 

where  n is  an  odd  integer.  Eq.  (21)  shows  that  the  closure  condition  (9)  is 
violated.  The  streamline  does  not  return  to  the  flat  "boundary  after  passing 
the  "bump.  It  is  obvious,  however,  that  the  streamline  could  be  made  to 
return  to  the  flat  boundary  by  the  introduction  of  discontinuities.  Since 
the  condition  (Al)  results  from  the  transient  motion  of  a sine  wave  boundary, 
the  general  configuration  is  symmetrical  about  the  mid-point  (X-^r  )•  The 

Mu 

discontinuities  must  bear  the  same  kind  of  symmetry;  and,  furthermore,  can 
only  occur  at  the  end  points  x = 0 and  x = TT . 

The  significance  of  these  seemingly  artificial  discontinuities 
probably  could  be  better  understood  from  a different  point  of  view.  Regard- 
ing (Al)  as  the  description  of  the  slope  along  a bump  in  steady  flow,  one 
sees  Immediately  that  the  shape  of  the  hypothetical  bump  must  be  a cosine 
wave.  If  an  up-and-down  displacement  of  the  cosine  wave  is  made,  two  dis- 
continuities are  brought  in  at  the  end  point.  These  discontinuities  being 
limited  to  a very  small  segment,  will  not  alter  the  condition  (Al) , but  will 
serve  to  change  the  streamline  directions  quite  abruptly  in  rising  from  and 
returning  back  to  the  flat  boundary. 
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To  ds tormina  the  off tot  of  these  discontinuities  on  tht  flow  by 


msans  of  the  linearised  theory,  a limiting  process  is  necessary.  Let  Bq. 


(Al)  first  he  modified  into 

<b  I •An  SLn  f*  Sin  nx  > x<"Tf 
Vq«o 

• 111  X— f - x **0 

S ( 'TT<k*rr+s 


- 0 


elsewhere 


where  £ is  a small  quantity.  The  magnitude,  An  . is  chosen  to  satisfy  the 


closure  condition  (9)* 


ft  , /TT 

/ ♦ I Ah  Sin  £ sin  n*  dx  * f dx*o 

ic  % J4  IT  > 


hencs 


” *? 


remembering  that  n is  odd.  Using  Xq.  (7).  the  increment  to  B(c()  becomes 


-fT  ♦ & 


**(«>'& ([(+ir  )(*£  *l**4*) 


/ .iu(n-ri)  mc\ 

7T«(''e  + e~  "e  ' 


Then  by  Xq.  (6) , we  have 


>d  I7  (x.)  - - 


Ah  r -ixiy  -tteS)l 

irrfttL  J 


Taking  the  limit,  5 "♦‘O 
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41  <b  (» 


Lit**  L 


L oc-  « 


-toe  IT 


'(<«)] 


a / z -tit-fr  \ 

a -*£L  (/•*•€  y 
ir* 

Saba tl tut log  Into  Xq.  (3)  • w*  find  that 


(A6) 


( Co*  * *♦  Co*  •tCTr-x)^*  -An  Seri  £»  fees  otic  * c*«e (fT-x)i 
TT  * n<rr  / * d<£ 


- /n  f^o[G(«)+cift)(^->«07 


(*7) 


whan  (A7)  la  added  to  Iq.  (19)  • tha  divergent  tarma  for  odd  h cancel.  The 
expreaaion  given  by  Bq.  (19*)  1*  thua  aatahllahad. 
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AFPBfDZZ  B 


Svaluation  of  Integrals  Appearing  in  the  Generalised  Force 


The  integrals  introduced  in  1%.  (38)  will  he  evaluated  in  this 


fir 

(1)  J$  J S,  (n*  ) Sin  m x clx 

fir 

(2)  J5 Ji(nim)a.j  5,(nx)c#s  mx  dx 

(3)  Tt>.»W  C,  (nx)S(f»  roxdx 

(4)  Jc  (n,m)^  J C,(nx)ct*  m x dx 

& 'l 

(5)  y^Stn  rk  in  d * 

fir 

(6)  / Sir>  rx  -^flk)dx 


(1)  Ts  (",™) 


By  the  definition  of  S,  (nx^we  have 


f " 


*4t 


^5 in  mxdx 
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Changing  the  order  of  integration, 

n 


PAGE  37 


r>?  x ax 


S*.n 


I f ’(  j r*i+ \ d."t 

mj  p-cosftvnr-— J+c^^^cosmTTj  — 


/—cos  rn *rr 
rw 


nmt  \ 

Sin  t cos  n Sir?  ± l 

* •t  j 


. r 


m 


| St  (mr)+£  Si  ( j S { (±  (n-w)'T^ j 


(Bl) 


where  the  upper  signs  in  the  last  term  are  used  when  n > m,  and  the  lower 


signs  are  used  when  n<  m.  In  particular, 

J^((«,o)a0 


(2) 

mx  d * 

after  changing  the  order  of  intergration.  The  integral  being  a divergent 
one,  a limiting  procedure  must  be  adopted.  Thus  for  n ^ m,  m ^ 0,  we  have 


In  a similar  manner, 

/v , (it-*)  /rn'rr)  . 

Jse  {f  +f  ^ 

'a  *•  'a  r o / 


dt 


r />< 


f nTT 

— if  cos  W TT  I 5J_rt  ~7T  t_  5i_Q_t 

* l -t 


1 


APiio.Fi  A^nr  ana  <;tpi  im  tprg  prqr apph 
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= “ “Vw~“ f CiOn''r'/‘fr)-tL  (<"♦•»>")  + 

flTe  f-Ci(  I'SIt)*  Ci  (‘  '**)&}] 


The  logarithmic  term  la  the  proper  limit  of  the  aum  In  the  curly  bracket, 
since  the  aaymptotlc  behavior  of  C£  la  known  (Ref.  2). 

Ci  (x) in  X , X *<  I 
^ being  Ruler’ a Conataat,  equal  to  1.78107 

Similarly  for  n - a,  n ^ 0 

/nir 

* nit 

* — )_+coS  tvrr  lio»  [ i -cot  el 

* n C-*o  J t 
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Using  the  asymptotic  "behavior  given  by  Eq,.  (B3) . we  finally  obtain 

»2L[-C£  (*/>*)+ fin  *'XniT]  (B4) 

Another  special  case  results  when  m * 0.  The  result  is  readily 
obtained  by  the  same  procedure, 


Si.  r»1T-|) 

(3)  Jct  («»»*») 

-If  n*  ) 

{(  -Jt  ^ 

•+  — 


ar-v ¥»'- 


xdx 


/n(iT-x)  -nk 

and  / 

by  their  proper  limiting  values.  The  limiting  value  of  the  combination  in 

the  curly  bracket  must  exist  because  it  Is  simply  f , though  it  is 

*nx 

artificially  split  into  two  terms.  The  integration  may  next  be  carried  out 
easily  after  a change  of  the  order  of  integration  as  was  done  for  J 
The  result  is,  for  njtn,  m^O, 
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for  r)  - m1  KM  ^ 0 


Ic  0vO  - - |c « (2 n ir)-2  Ci  (n it) ft 


(B7) 


and,  for  KM  - O 

Jc  (Nn,  o')  - O 

wTc  (»,"0 


7 Avn)  involve  terms  which  appear  in  T . The  result  is, 
ez  ' J ^>1 

^«r  m \ 0 

3c^(h,  noV  ~C^^l>t7T  Si (±(n-rri)rr  )| 


(B8) 


where  in  the  last  term  the  tipper  signs  are  used  when  n > a and  lower  signs  are 
used  when  a<,  m.  Tor  m * 0,  J,  (n  6^-0  • 

/ r 

(5)  J sir*  r*  JL-  d* 

r« 

dx  - J 5tn  rx  x dx  — 

j Sin  rx  JL  (MT- *)  dx  = 


fir 

J Sin  rx  Jin  x 


ir-x 


, \ /V 

\l  4-  Cos  riT)  J ?in  rx  An 


dx 


(B9) 


which  result  is  also  obvious  from  symmetry  considerations.  Integrating  by 
parts, 


AEROELASTIC  AND  STRUCTURES  RESEARCH 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY  . 
Department  of  Aeronautical  Engineering 


CONTRACT  NO.  15  ori -07833 

sa  06*1-259 


PAGE  41 


P /If 


Sin  ry  2.w  x Jx  * — “ * j + jr  f 

• /• 


--/rfccs  r Tf  Jtn  Tf  - C+( rir)-  ( r£)f  ] 


«***• 


(BIO) 


Bence • 


J Sin  r*  1.X  J*  = — I *-  f.l?  po  r«  ILn-n -CiOir'ii-Jl^yf] 


l + C.t 


r-TT  Y-C!(rtT)7 


(m) 


Si"  r*  Cx'^Ta) 


igdfi  from  wye  try  conalderetlone , 


jfs>»  ^(4-Tr-x)Jx=(,+c#t 

(i  + cnsrTrjSilr^) 


(B12) 
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Ic  fci(zmr)-2Ci(nir)t  ■ (b?) 

and,  for  W ' O 

JCi  Qrr,  <0  = o 

wl  Oo,*0 

^2 

j:  <v)  involve  terms  which  appear  in  T . The  result  is, 
cx  j 

*r?  \ 0 

|-c2°^Tr  jsi(<"«-’0'>r)?  S;(±(r,-m»j  <M) 

where  in  the  last  term  the  upper  signs  are  us&d  when  n > a and  lover  signs  f re 
used  when  n<;  m.  Tor  a * 0,  (n,6}*0  • 

sV 

(5)  j sir*  rK  &r>  J<__  cjx 
TT 

/Si  n rx  JLr»  — dx  = I ^in  rx  * dx  - 

TT-X  4 

f Si  r>  rx  hr>  0 Tf - X ) d X = 
i \ 

o 

( f +•  Cos  r rr^  sinrx  in  x dx 


which  result  is  also  obvious  from  symmetry  considerations.  Integrating  by 
parts , 
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f7*  Sin  r.-i*  x Jx  *=  *• 


rx  fin  x I . i f 

r rr 

• Vo 


s — - fr[coS  — 6»  flry>rJ  — { </•*  £ ~ J 


f-**® 


(BIO) 


Bence, 


'Tpr  — 

y S*r>  rx  JLnX  Jv  = — JJ^Cr^rnT  r „ jin  -fT  - Ci  ( r tr)  -f-J0a  fr] 


r 'TT  Y ' C i (r  1T)J 


T 

/ S(n  f*  (*-£*)  dx 

v# 


Again  from  ayn  try  eonaidemtlo&s , 


/’sin  rx  (^-^._x')<lx=(/+CoS  rir)pM 

(i  + c«p,n')5i(fff) 


2-tiS-d; 
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